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Abstract 

It is pointed out that the average semi-inclusive particle phase- 
space density at freeze-out can be determined from the coincidence 
probability of the events observed in multiparticle production. The 
method of measurement is described and its accuracy examined. 

1. Recently, several methods were proposed which allow to estimate the 
average of the single-particle inclusive phase-space density produced in ultra- 
relativistic heavy ion collisions PP-jS]. This quantity is useful in discussions 
of the equilibrated systems and therefore such measurements open possibil- 
ities to verify the expected presence of the thermal phase-space distribution 
at freeze-out [U] and/or search for more exotic phenomena jjj- 

In the present paper we propose an extension of these studies by in- 
cluding, in addition, the exclusive and semi-inclusive M-particle phase-space 
densities. We show that their averages can be estimated from measured co- 
incidence probabilities of the multiparticle events observed in high-energy 
collisions. The information one may gain from this approach is complemen- 
tary to that obtained from single-particle inclusive measurements. In partic- 
ular, it gives an insight into the correlation structure of the final state of the 
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collision (including both dynamic and Bose-Einstein correlations), a feature 
which is ignored in the single-particle inclusive measurements described in 
PQ-[5]. Furthermore, one can show [8| that the average semi-inclusive phase- 
space densities are closely related to the second Renyi entropy and thus 
their measurement allows to estimate a lower limit of the true (Shannon) 
entropy of the system (without assuming the thermodynamic equilibrium). 
Needless to say, a comparison of the measured average semi-inclusive phase- 
space densities with expectations for the thermalized systems would be very 
interesting indeed. 

As emphasized in [UIH], the phase-space distribution of particles produced 
in high-energy collisions is not a precisely defined quantity. Apart from 
the standard problems with the uncertainty principle, one has to take into 
account that particles may be produced at different times JI|. Following 
Bertsch P, we shall ignore this problem and assume that all particles are 
created at the same time. That is to say, we are considering a time-average 
of the density [TJ |3| . 

To define the average semi-inclusive phase-space density, consider a collec- 
tion of events in which exactly M particles were observed in a given region of 
the momentum space. We shall call them M-particle events (independently 
of how many particles were actually produced) 1 . These events can be de- 
scribed by the normalized M particle phase-space distribution Wm{X,K), 
with X = Xi, Z M , K = K^-\ K( m >. The corresponding particle phase- 
space density is D M (X, K) = MWm(X, K) and thus 



gives the average phase-space density of the M particle system. 

It should be emphasized that, as is clear from this discussion, the phase- 
space density D M , describing the semi-exclusive distribution, refers only to 
particles actually measured in a given experiment and in a given mometum 
region. It gives no direct information about the particles which are not regis- 
tered in the detector. To obtain information on the phase-space density of all 

1 This is terminology often used in experimental description of multiparticle processes. 
The proper technical term is the exclusive distribution if all particles are observed, and 
semi-inclusive distribution if besides a given number of observed particles there is an 
unspecified number of other particles. This should not be confused with inclusive M- 
particle distributions. 
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produced particles additional assumptions (e.g. thermodynamic equilibrium) 
are necessary. 

Note that the average semi-inclusive phase-space density averaged over 
all particle multiplicities is simply obtained from (??), using 

« D »= J2 p (M) < D M > /77c (2) 

M 

where P(M) is the multiplicity distribution. Therefore from now on, to 
simplify the discussion, we shall only consider the case of a fixed multiplicity. 

Our method is based on the observation that, for a rather wide class 
of models of particle production, the quantitity C, defined in (??) 2 , can be 
approximated by the measured coincidence probability C exp of the events 
with M particles, defined as jTOJ [TT] 

where N2 is the number of the observed pairs of identical events and N is the 
total number of events. N(N — l)/2 is the total number of pairs of events 3 . 

It is clear that, since the observed events are described by the particle 
momenta which are continuous variables, Eq. (??) is not directly applicable: 
a discretization is necessary. Then one can define the identical events as those 
which have the same population of the predefined bins and thus counting 
of coincidences becomes straightforward 4 . The counting of identical events 
obviously depends on the binning, so the procedure is ambiguous [TOJ HU H31 
ITo] . In order to obtain a viable estimate of the average particle density, we 
thus have to select the binning in such a way that the result of (??) is as close 
as possible to the exact value of C which, as seen from (??), gives directly 
the particle phase-space density. 

In the present paper we argue that for a fairly large class of physically 
sensible models, one can find an adequate binning procedure and thus to 
determine rather precisely < D > by measuring C exp , i.e., by counting the 
number of pairs of identical events. The method turns out to be particularly 
suitable for large systems and thus may be useful in heavy ion collisions. 

2 To simplify the formulae we shall from now on omit the index M in all quantities. 
Since we are discussing solely M-particle events, this should not lead to any confusion. 
3 Formula (??) was first suggested, in a different context, by Ma |12j . 
4 A detailed description of this procedure was given in P2J and applied in [T3] , 
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In the next section the discretization procedure is described in some detail 
and the corresponding formulae for C exp are written down. The phase-space 
density < D > and its relation to C exp are discussed in Section 3. Our 
conclusions and outlook are given in the last section. 

2. In this section we discuss how the discretization procedure afects the 
definition (??) of the coincidence probability. To this end we first express 
the C exp given by (??) in terms of the momentum distribution of M particles 
w(K)=w(KW,...,KM). 

Consider a set of discretized events constructed by dividing the particle 
momentum space into J rectangular bins of volume 



Uj = (A x A y A z )j ; j = 1, J.ma (4) 



Then the probability to find a particle in bin u^, another one in bin Uj 2 , etc. 
is 

M 

P(ji,32,..,jM) = II u jm <w(K$;\...,KW)>l8l (5) 



m=l 

where 



w (K W K {M) ) 
w V n ' 3m ) 

= n(^J -1 / dK^...( dK^w(K^,...,K^).ls2 (6) 

m=l ^3M 



<w(k):',...,ky;:>) >-. 

M 



Note that the bins cu^, ■■■,Wj M do not have to be different. 

Thus the coincidence probability as measured by the formula (??) is 



jl 3M 

M 

= E-"E m*iJ 2 [<™(4V.,K£) >]^3 (7) 

jl 3M m=l 

The first equality follows from the observation that sampling a series of events 
is the Bernoulli process and thus probability to find, after iV trials, ni, nj 
events in configurations {1}, ...{J} is 

B( ni , ...nj) = -J^—(p i yK..(P J ) n n?>2 (8) 
n 1 \...n J \ 
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From this formula it is not difficult to see that 

<iV 2 >. £ £^^B {nu ...n J ) = ?Wf»£(P i ri 3 3 (9) 

ni,...,njj=l A / j=l 

The question now is: how to select the bins u) 3 - to obtain a result as close 
as possible to C giving the average value of the particle phase-space density 
< D > [c.f. (??)]. This is discussed in the next section. 

3. To analyze the relation between C exp and C we consider the M-particle 
phase-space distribution of the general form 

W(X, K) = 1 — G[X/L]w(K)ma (10) 

with X/L = (Xi — Xi)/L x , (Z M — Z M )/ L z , K = Ki, K M . The function 
G satisfies the normalization conditions 

J d 3M uG(u) = 1 -»• J dXG(X/L) = (L x L y L z ) M ; 
Jd 3M u Ui G(u) = -> <X u Y u Z l >=X l ,Y u Z l ; 
J d 3M u( Ui ) 2 G(u) = 1 -> < (Xi -Xi) 2 ,... >= L 2 X ,...W (11) 

The first condition insures that w(K) is the observed (multidimensional) 
momentum distribution 5 , the second defines the central values of the par- 
ticle distribution in configuration space and the third defines L x ,L y ,L z as 
giving root mean square sizes of the distribution in configuration space. Both 
sizes and central positions may depend on the particle momenta 6 . The form 
of the function G describes the shape of the multiparticle distribution in 
configuration space. 

Ansatz (??) for the time-averaged phase space density is satisfied in a 
large variety of models 16 . 

Using (??) we obtain from (??) 

C= (2*)*" / d™K[ W { Kl ,...,K M )T [ ^|[G(I/I)f = 

= (2ng) 3M [ d 3M K[w(K 1} K M )f - \ 138 (12) 

J \ L >x L >y L >z) 

5 By definition, w{K) = J dXW{X,K). 

6 They may be also different for different kinds of particles. 
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with 

g 3M = J d 3M^ G ^2 l39 (13) 

The constant g depends on the shape of particle distribution in config- 
uration space. This dependence is, however, rather mild. For example, we 
obtain g~ l = for Gaussians and g~ l = 2y / 3 for a rectangular box. 

In the discretized form, (??) can be written as 



M < , 

C = (2ng)™ £ n , r Tl \ < Kwm] 2 > MO (14) 

h,...,j M m=l \^x^y^z)j m 

where 



M . , 9 

IlKJ- 1 / dK^... dK^{w(K^,...,K^)} 2 l36 (15) 

m=l w Ji J ^Af 



and (L x L y L z )j m is a suitable average of (L x L y L z ) over bin j m . 

Comparing (??) with (??) one sees that to have C^ p as close as possible 
to C, the best volume of the bins is 

u, jm = (A x A y A z ) 3m = {2n f Ml (16) 

\- L >x- L >y- L >z)j m 

One sees from this formula that depends crucially on the volume of the 
system in configuration space. One sees, furthermore, that with this choice 
of uj the coincidence probability C exp , determined by counting the number 
of pairs of identical events [c.f. (??)] 4 , is related to C, giving directly the 
average particle phase-space density < D > [c.f. (??) and (??)], by the 
formula 



Q _ £<exp '52ji,---,jM ^ [ W ji,---,jA'l] ^ ^jo (17) 
SjlvjAfi^ W 31,-,3M > ] 2 

It is thus clear that the accuracy of determination of C increases with 
increasing volume of the system. Indeed, for a volume large enough, the bins 
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defined by (??) are small and the ratio on the R.H.S. of (??) approaches 
unity. For smaller volumes the method is less accurate but one may try to 
estimate the correcting ratio from the (measured) single particle distribution. 
4. Several comments are in order. 

(i) One sees from (??) that the optimal size of the bin does not depend 
on the average position of the particles at freezeout. This implies that the 
momentum-position correlations induced by the i^-dependence of X do not 
influence significantly the measurement of the coincidence probability. 

(ii) It is also seen from (??) that only the volume of the bin ujj m = 
(A x A y A z )j m , but not its shape, matters in the determination of the optimal 
discretization. One can use this freedom to improve the accuracy of the 
measurement by taking bins large in the directions with weak momentum 
dependence and small in the direction where the momentum dependence is 
significant. 

(iii) One may improve the accuracy of the measurement by estimating 
the ratio {Ej u ...,j u < [ w h,...,j M } 2 >}/{Ej lt ... 

j'jvit^ w ji,---jM }• This may 
be possible if the momentum distribution of particles is measured with good 
accuracy. 

(iv) Our analysis can be applied to any part of the momentum space. 
This allows to measure the local particle density in momentum space, av- 
eraged over all configuration space. In case of strong momentum-position 
correlations, the selection of a given momentum region can induce, however, 
a selection of a corresponding region in configuration space. 

(v) The accuracy of the measurement depends crucially on the correct 
estimate of the size of the system. Information from HBT measurements 
should allow to determine the parameters L x ,L y , L z and -at least in principle- 
also the shape 7 of the function G(u) (some procedures are described in [21 El)- 
Therefore good HBT data are essential for a successful application of the 
method. 

(vi) The presented analysis of the discretization procedure can be gener- 
alized to higher order coincidence probabilities This opens the way to 
a determination of higher Renyi entropies [H] and then, by extrapolation, to 
obtain information on the Shannon entropy of the system [TTj . 

In conclusion, we propose to estimate the phase-space density of particles 
produced in high-energy collisions by measuring the coincidence probability 

7 Needed to evaluate the parameter g [cf. (??)]. Fortunately, as we already noted, the 
sensitivity of g to the shape of G(u) is rather mild. 
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of the observed events. The accuracy of the determination of the coincidence 
probability by counting the number of the identical events fUJ E] was anal- 
ysed for a large class of physically sensible models. It was shown that the 
accuracy improves with increasing volume of the system and, therefore, the 
method is particularly suitable for heavy ion collisions. A formula giving 
the optimal discretization method in terms of the size of the system in the 
configuration space [Eq. (??)] was derived. 
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